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The consequence of dynamical mass generation on the singularity structure of prop- 
agators is discussed. First the phenomena of dynamical mass generation is discussed in 
the framework of Euclidean gap equations, then a possible Minkowski solution is looking 
for. The examples are reviewed and studied for several models: Yukawa, QED, QCD and 
Wess-Zumino. It is argued that the absence of propagator pole goes hand by hand with 
the nontrivial solution for mass function. The consequences are discussed. 
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■ 1 Introduction 

In perturbation theory the S-matrix elements are composed from Greens func- 
tions which are calculated and regularized/renormalized order by order according 
to the standard rules. By construction, the correct analyticity of the S-matrix is 
ensured by the correct analyticity of Greens functions. As a consequence of as- 
sumed analyticity various Greens functions satisfy known integral representations. 
O |' As a simple but important example recall the spectral representation of fermion 

propagator 
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and the appropriate dispersion relation for selfenergy function: 

= f d MH + PsH (2) 
j p-* - uj + it 

where two (in parity conserving theory) cr's are called spectral functions. These 
are uniquely determined by knowledge of selfenergy absorptive parts p's (note 
S^ 1 (j>) = i>— m— £(p), m is a fermion mass and we do not address the question of 
renormalization here) . In the case of small coupling constant the later functions can 
be read from the perturbation theory series of Feynman graphs. In many practical 
cases the use of perturbation theory is not only very useful, but almost unavoidable 
tool for meaningful predictions in particle physics. 

In this paper we explore the analytical properties of propagators in strong cou- 
pling theories where nonperturbative approaches are necessary and we deal with 
the quantum field models with dynamical mass generation (DMG) . For this purpose 
we consider no mass term of given field in the Lagrangian quantum field theory, 
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but assume nontrivial solution of the mass function. In the case of fermion field, 
the mass function M is conventionally defined as 



where F(p) is usually called fermion renormalization function. The identification 
with the previous definitions is straightforward. For a review of strong in- 

teracting models with the effect of DMG studied in the context of gap equations 
see p. 

In this paper we demonstrate on few model examples that the correlation func- 
tion of the field with DMG does not posses a pole on the real positive semi-axis 
of complex momenta. In other words the 'would be particle' associated with the 
field cannot live on-shell. As an hopeful convincing nonperturbative tool we use the 
formalism of Schwinger-Dyson (gap) equations and/or certain numerical analysis 
of some recent lattice QCD data. Recall, the gap equations are exact consequences 
of quantum theory and they are derivable from the appropriate path integral [2] ■ 

The theories with phenomenon of DMG play an important role in the particle 
physics. First, let us mention the Standard Model. There, the dynamical chiral 
symmetry breaking is clearly inherent (although approximate) property of light 
quark sector in QCD, while the tree or 'bare' masses are generated via usual Higgs 
mechanism. Furthermore, as a very tempting alternatives to the SM the Technicol- 
ors (or extended Technicolors) models were developed (for a review see [3 ) . These 
typically Higgsless models do not require non-zero expectation value of scalar field 
and the electroweak symmetry SU(2) w *U(1)y /U(l) em is broken dynamically. This 
is associated with formation of fermion condensates, which is typically driven by 
a new, possibly non-Abelian strong force and hence it reminds global chiral sym- 
metry breaking SU(2)l * SU(2)r/SU(2)v in u,d quark's QCD sector. Having no 
direct experimental evidence for Technicolor scenario the masses of a techniparticles 
should be scaled up to (at least) several TeV. For a possible indirect experimental 
evidence and some recent constrains on technifermion models see the papers 0], 
. Let us note for completeness here that the role of scalar field and the usual (but 
now strong) Yukawa interaction can be reconsidered for this purpose jEJ,[I|- 

Further, there is another set of interesting models which suggests the radiative 
origin of (at least the lightest) massive fermions. The so called Zee model [E],[II] 
assume the existence of an additional scalar SU{2) doublet with the ordinary La- 
grangian mass term. The origin of small neutrino masses should be economically 
explained as a pure radiative corrections effect that follows from the Yukawa cou- 
pling constant matrix. Regarding contemporary measurements of neutrino mixing 
the Zee model is still one of attractive models for neutrino mass matrix. In order to 
enhance the reliability of the model at high energy scale, the Zee model has been 
embedded into the supersymmetry |1(J| . 
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2 Using Standard Methods 

On an example of the simple model we review or mention the standard methods 
usually used in the literature. For this purpose we choose the model already investi- 
gated by Bicudo in 1 1 1| . This is in fact simplification of Schwinger-Dyson equation 
for fermion in gauge theory and/or in the Yukawa model with the (nonhcrmitian) 
interaction igipipcj) where only the scalar field mass term m 2 (j> 2 /2 is explicitly in- 
cluded in the Lagrangian. For further simplification the ultraviolet divergences are 
not addressed here and the interaction kernel of gap equation is regularized by 
a Pauli-Villars term. Having the integral kernel finite, we approximate the field 
renormalization function, F = 1 in the equations. Then the 'mass gap equation' is 
a non-linear integral equation for the mass function M, 



d 4 q_ 
? 



T(q,l)G A (p-q)S(q), 



Ga(x) = o - — 9 .9 ; ■ (4) 

where T is a boson-fermion interacting vertex, G A is boson propagator with 
Pauli-Villars regulator. The trace is taken over the dirac indices. Boson propagator 
as well as the vertex T satisfy their own Schwinger-Dyson equations. For the further 
technical simplification we take them in their bare form, i.e T = 1, G _1 = p 2 — m?. 

Chiral symmetry breaking is associated with a nontrivial solution for a mass 
function M. Depending on the details of the model it typically appears above a 
certain critical value of the coupling a c , where as usual we denote a = g 2 /(Air). 
Hence to solve the problem of dynamical mass generation is nontrivial task since the 
kernel of eq. involves (potentially complex) singularities while the value of the 
coupling, a c ~ f suffer from the usage of known standard perturbative technique. 
Recall here that a c — 7r/4 m Euclidean formulation of presented model when 
mt = and when a hard cutoff regularization is used. 

Substituting the Ansatz JIJ into eq. J3J we get: 

Mb 2 ) = e!^-Jd^ 2 a2iW ]. G A (p- q ) (5) 
J (2tt) 4 J p z - lj + le 

In the standard manner we can obtain the resulting dispersion relation for the mass 
function M 

M(p 2 ) = (to / sM . ' 
y ' J p 2 -uj + te 

a 



= J^j J dxa 2 {x) [X (u;ml,x) - X (u] A 2 ,x)] , 
■o(«;a,6) = ^Se^^^f), ( 6 ) 
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where X(lu, a, b) is well-known triangle function 

X(u), a,b) = (u - a- b) 2 - Ayz = u 2 + a 2 + b 2 - 2uja - 2ub - 2ab 



The imaginary part of propagator and the absorptive part of selfenergy M are 
related through the following Unitary equation 13 , 



Ps(P 2 ) , C 2 v {p 2 ) 



1 I 1 
p (T 2 [P 



c s {p 2 ) c 2 s (p 2 ) 
c s {p 2 ) = p 2 + ReM 2 ( P 2 ) + 7r 2 p 2 s (p 2 ) 
c v (p 2 ) = -2ReM(p 2 ). 



(7) 



which, as can be deduced from the Heaviside step functions in eq. 10, are nonzero 
above the threshold. In eq. J7J the above index c means the continuous part of 
the spectral function. From its derivation it follows that the eq. J7J) requires the 
principal value integration. This can be performed in analytical closed form: 



ReM(p 2 ) = P. J = / dx ^(x) [j(p 2 ,a;,TO 2 ) - J(p 2 , x, A 2 )] , 



p* — LO 
where we have labeled 

©(-A P )v^ 
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(9) 



where we have introduced the following abbreviations: 

X P = X(p 2 ,y,z), 

\/% = \v-z\, 



(10) 



Note that having m& different from zero one can exclude the possibility that the 
spectral function a is purely continuous. The inverse of propagator becomes com- 
plex from the threshold (T + mj) 2 while T is a lower boundary of the spectral 
integral ijT)). i.e the point where the propagator becomes complex. This is clear 
contradiction and the Unitary Equations have no solution in this case. Hence we 
assume usual following form 

<7 2 (x) = rS(x - to 2 ) + 0-2(2) (11) 
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where r is a propagator residuum and m is a pole mass. These two propagator 
characteristics can be selfconsistently determined as follows 

1 



m z — x 

As the second and in fact most spread technical method of solution we make 
our calculation in the Euclidean space. In the Euclidean formalism the analyticity 
domain of propagator functions should be the same as in the previous spectral- 
Minkowski treatment. Assuming no poles in the first and third quadrants of complex 
momenta one can deform the integration contour and perform the so-called Wick 
rotation iq — > q 4 , where 54 is the fourth component of Euclidean momenta q E = 
g 2 (and similarly for p). Choosing then convenient angular parametrization 

qi = qE sin sin sin 77 

<Z2 = 1e cos (f> sin 6 sin 77 

<Z3 = Qe cos (9 sin 77 

<74 = tecosry (13) 

such that (j>e — Qe) 2 — P% — %Pe<1e cos 77 + q E one can integrate over all angles 
which leave us with only one dimensional integral equation for the function M 



M (v 2 ) ~ r da 2 M(a 2 ) ^ + - + ^ X{p ^ q ^ < ] = ^ X(p *> q ^ A " ] (U) 

(see also [H]). The equation (|14|l can be easily solved numerically by the method 
of iteration. 

Assuming uniqueness of the spectral representation for fermion propagator and 
using the theorem of uniqueness of the analytical continuation we should find that 
the function M(q 2 ) obtained as a solution of the Unitary Equations (JJJ must be also 
the solution of the equation (|14|l for M(q E ) at spacelike momenta. The disagreement 
implies the 'unexpected' position of Greens functions singularity. We should stress 
here that the spectral and the Euclidean solutions actually coincide in weak coupling 
quantum field models with explicit mass presented [13] ■ 

Before the presentation of the results let us mentioned further possible method 
of solving gap equations in Minkowski space. This is based on a partial neglecting 
of the time momentum component in the kernel of the equation. Typically only 
the assumed dominant part of propagator poles is considered. These and similar 
approaches were mostly used in the calculation of bound state equation, while for 
a given model some review can be found in However, the dynamical chiral 
symmetry breaking is deeply nonperturbative -strong coupling- phenomenon and 
we believe it is difficult to have such an approximations under safe control and 
we do not employ this method. In what follows the numerical iteration of the 



Czech. J. Phys. -1949 (1) 



X 



5 



V. Sauli 




1(T 10° 1(T 10 4 10' 10" 10' ! 10° l(f I0 4 



2 2 
"P "P 



Fig. 1. Dynamical mass as obtained in the Euclidean (E) formalism and spectral (S) 
formalism. Left (right) panel shows the result in log-log (log-linear) axis scale. Each line 
is labeled by the coupling constant. 

equations offers stabfe convergent and accurate solutions in both approaches. We 
take the boson mass mj, = 1 and all the results are scaled in these units in our 
actual calculations. The square of the cutoff is A 2 = 50. In the Euclidean theory- 
such constellation leads to the expected dynamical chiral symmetry for the coupling 
larger then a c ~ 6. The infrared value of generated mass, say M(0), has a drastic 
dependence on the coupling constant. We did not attempted to parametrize this 
behaviour but recall here the paper |12| where, within a zero boson mass the authors 
of |12| have found the well-known exponential behaviour (this so called Miransky 
scaling is explicitly shown in the next section). The appropriate results for the 
function M(p 2 E ) are displayed in fig.l. In what concerns the behaviour of the mass 
at ultraviolet, Bicudo shows that 

MCr*)- ^ 2 -"^^ 2 , (15) 

where the fermion condensate is defined as usual < ipip >=gTr S. This power 
behaviour seems to be confirmed by our numerical analyzes. 

Unfortunately, this is not what we do observe from the solution of the unitary 
equations. First of all let us mentioned that there is no nontrivial critical coupling. 
Stable solution for the functions a c (p s ), residuum r and pole mass m is obtained for 
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any non-zero coupling constant a. The mass function for spacelike p 2 is then calcu- 
lated through the appropriate dispersion relation and compared with M(p 2 E ). These 
results are added to the fig. 1 for several values of coupling constant. From this we 
can see that these two solutions do not agree for the same coupling. Furthermore, 
the scaling M(0) is approximately linear for our 'would be' Minkowski solution. 
Also the ultraviolet behaviour differs, the 'spectral solution' for M is driven by the 
first power of the inverse square of momentum, i.e. M(p 2 ) — > const/p 2 . Note, the 
possibility that this is also a solution of ( non-unique for a while) the gap equation 
i|14|) has been excluded by simple substitution and comparison. 

Using analyticity arguments above we can conclude that the structure of corre- 
lation functions at the complex momentum plane prohibits naive backward Wick 
rotation ip^ — > po and the propagators does not posses naively expected spec- 
tral (Lehmann) representation. Here we stop our discussion and leave the correct 
Minkowski treatment for a future interesting investigation. 



Considering explicitly massive electron, the analytical structure of strong cou- 
pling QED electron propagator was subject of the initial study It was found 
that above the certain value of the coupling constant the branch point of the fermion 
propagator disappears from the real timelike semi-axis. This was also confirmed by 
the later author's study made in the paper [T5| . Furthermore, based on the solution 
of Unitary Equations presented in the previous section, it was argued that spectral 
representation Q is absent in these circumstances (for a further review of various 
Euclidean solutions in QED and for a review of similar spectral techniques see the 
list of references in ^3]>GH])- It is very interesting fact that the above mentioned 
value of the coupling is exactly identical with the critical coupling a c = ^ which is 
true at least in the ladder approximation of Schwinger-Dyson equation for S. 

Remind at this place the trivial fact that the existence of free particle moving is 
necessary associated with a presence of propagator pole at some momentum square 
p 2 = m 2 > 0. Hence, if no fermion mass term is considered in the Lagrangian 
and simultaneously free massive fcrmions are expected then the naive perturbative 



dispersion law E = \ffp must be nontrivially transformed into E = ^Jp 2 + m 2 . 



Here the 'gap' is naturally identified with the pole mass m p , which must be a 
solution of our gap equation. On the other hand the absence of a real propagator 
pole could be a hint for a confinement or resonant behaviour. 

From the above arguments one can conjecture that the confinement is somehow 
intimately related with the dynamical chiral symmetry breaking. At the introduc- 
tory section we enumerate several models that propose the dynamical origin of the 
particles in the nature. Clearly, having the leptons unconfined it is desirable to look 
for an interaction which spawn DMG but avoid a confinement. In a gauge theory, 
the behaviour of the running coupling associated with the proposed new strong 
interaction plays the key role. As an simple example let us consider massless gauge 
theory where the coupling is not running but is approximated by a constant value. 
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For further simplification the fermion-antifermion-photon vertex is approximated 
by the bare vertex 7^. In this so called Ladder approximation the fermion gap 
equation reads: 

2 f d 4 k 



Z(p)=£-ie*J ^*Gr(k-ph»Stoh v (16) 

where is a free massless gauge boson propagator in Landau gauge. 

Recall again the most conventional strategy: after the Wick rotation and the 
angular integration the gap equation in Euclidean space is formulated: 

, 3a f A ' „ 2 2kl Mix) 

M(p E )=m + — dk% — — = - 2 ■ -y- - (17) 



47r7 * p 2 E + k 2 E + ^I—pT k 2 E+M 2 {kE) 

where m = when studying chiral symmetry breaking and where we have intro- 
duced an Euclidean cutoff regulator A as the upper boundary of Euclidean integral. 
Note that now, contrary to previous case, F — 1 as the exact solution because of 
ladder approximation. 

Like in the previous section the second alternative possibility is to assume spec- 
tral representation JIJ and to solve the Unitary Equation directly in Minkowski 
space. The evaluation of selfenergy and its absorptive part p s is straightforward 
(for details see for instance ^3]- Note only that the eq. (JJJ is exact now since 
p v = in ladder approximation. 

As a first step, mainly due to the pedagogical reasons, we exhibit the agreement 
of both mentioned approaches and present solutions for small (in fact subcritical) 
couplings in an explicit chiral symmetry breaking, i.e. for non-zero to . For this 
purpose we prefer to show the renormalized results obtained by the manner firstly 
used in the paper [TJ]. The 'spectral' solutions at the time-like regime are displayed 
in Fig. la. The maximum cusps in the real parts of dynamical mass function M(p) 
correspond to the thresholds, i.e. to the physical mass of the electron. The contin- 
uous imaginary part of M(p), i.e. iTp s (p 2 ) is the result of the solution of Unitary 
equations. These are displayed in Fig. 1 (left panel). The comparison of Euclidean 
and spectral solutions is shown at the right panel of Fig. 1. Solid line curves repre- 
sent the Euclidean solutions, the dot-dashed lines are the mass functions obtained 
via spectral treatment. The coupling constants a = e 2 /47r decrease from top to 
bottom. For completeness note that zero momentum subtraction renormalization 
scheme was used. At this place we should stress that only the pole mass, defined 
as Trip — Yi(m p ), is the observable ( it should be independent on the gauge and 
on the renormalization scheme) . In addition we switch off the electron mass in the 
Lagrangian. In the Euclidean formalism the properties of DMG are well known in 
this model. The mass function in the infrared obeys the Miransky scaling |17j 



M(0) ~ Ae^V 11 ^ (18) 
for a > a c — ?, while for a < a c chiral symmetry is preserved and S = K exactly. 
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Fig. 2. Fermion mass functions in the case of explicit chial symmetry breaking at timelike 
(left panel) and spacelike (right panel) momenta as described in the text. 



Dynamically generated mass is displayed in the figure 2 (solid lines). The ap- 
propriate curves are labeled by the coupling constant a. The results for dynamical 
mass M in unquenched QED are added for interest to. The polarization effect were 
taken selfconsistently into account and we have used the bare vertex in this case. 
The Schwinger-Dyson equations in unquenched QED have been solved in several 
papers |181 1131 1191 I16| . The results presented here follows the numerical treat- 
ment described in T3J. What is the timelike behavior? It was already shown by 
Fukuda&Kugo in 1976 14 that there is no evidence for the propagator pole in 
ladder QED. The function p 2 — M(p 2 ) never cuts the zero axis for a > a c . This 
observation is valid even for explicit chiral symmetry breaking case and it has been 
confirmed in the paper |15| again. Yet more, no evidence for solution of our unitary 
equations at supercritical regime a > a c was found (the same is observed in un- 
quenched case). It is probably strong hint (but of course not a proof) for a changes 
of propagator's analyticity domain and we argue that fermion propagator does not 
posses usual spectral representation in Ladder QED. We do not investigate the 
analytical structure of the fermion propagator in ladder QED further. Recall, that 
one can expect some but not radical influence from vertex corrections which was 
neglected in our study. 
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Fig. 3. Dynamical masses in ladder QED (solid lines) and beyond. 

4 QCD 

The Quantum Chromodynamics is the only experimentally studied strongly 
interacting relativistic quantum field theory. This non-Abelian gauge theory with 
a gauge group SU (3) has many interesting properties. The dynamical spontaneous 
breaking of chiral symmetry explains why the pions are light, identifying them with 
the pseudo-Goldstone bosons associated with the symmetry breaking of the group 
SU(2) L x SU(2) R to SU(2) V (in flavor space). Asymptotic freedom [201121] implies 
that the coupling constant of the strong interaction decreases in the ultraviolet 
region. For less than 33/2 quark flavors the QCD at high energy becomes predictable 
by the perturbation theory. However, in the infrared region the Perturbation theory 
does not work and non-perturbative techniques have to be applied. 

One of the most straightforward non-perturbative approaches is a solution of 
the gap equations for QCD. The extensive studies were undertaken for a quark 
gap equation, based on various model assumptions for a gluon propagator. These 
approximate solutions, often accompanied by a solution of the fermion-antifermion 
Bethe-Salpeter for meson states, have become an efficient tool for studies of many 
non-perturbative problems, e.g., the chiral symmetry breaking, low energy elec- 
troweak hadron form factors, strong form factors of exclusive processes, etc (see 
reviews [221 EH El EEj and also recent papers ffi\ 12*5] ) . 

However, to take gluons into account consistently is much more difficult than 
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to solve the quark SDE alone. At this place we briefly review the symmetry pre- 
serving gauge invariant solution obtained by Cornwall two decades ago |29| . To our 
knowledge this is the best published example, in which the behavior of the QCD 
Green's functions in the whole range of Minkowski formalism is addressed within 
the framework of the gap equations. 

Cornwall reveal the resulting gluon propagator which is gauge invariant by the 
construction and the gluon form factor does not depend on the gauge fixing pa- 
rameter, although it was used at the beginning of the author's calculation. In his 
approach the Green's functions satisfy the naive Ward identities which is a direct 
consequence of gauge and pinch technique introduction. For the details see the 
original paper |29| . It was shown that the resulting gluon propagator can be rather 
accurately fitted in spacelike regime as 

q»q v 



G^{q) A B = S AB G(q) \-gV + 
G~\q) = [-q 2 + m 2 (q 2 )]bg 2 ln 

"7 /-g 2 +4m 2 \ 
ln { ^ ) 



q 2 



q 2 + 4m 2 (q 2 ) 



It 2\ 2 

m (q ) = m 



A 2 
12/11 



where g is a gauge coupling, b is one loop rcnormgroup coefficient, m ~ [1.5 — 2]A, 
A has the meaning of Aqcd- The Latin (Greek) letters stand for a color (Lorentz) 
indices. One should mention here that this propagator was successfully used in some 
recent phenomenological calculations of r decay, meson form factors, etc.. , see for 
instance [501 15T1 152] . 

In recent papers |2I3 1221 EU ESI EH studies of the coupled SDEs for gluon 
and ghost propagators in the Landau gauge in Euclidean space were performed in 
various approximations. However the non-Abelian character of the QCD makes it 
difficult to convert the momentum SDEs into equations for spectral functions. We 
were not able to do this yet since the Fadeev-Popov ghosts have to be included 
|57] in a class of Lorentz gauges. The main obstacle is the ghost SR, mainly due to 
singular behavior at zero momenta (see for instance ^26 ). 

Hence, instead of solving the SDEs the following generalized spectral represen- 
tation 

00 

G^)=[^-p^, (20) 
J q z - lo + ie 


was used to fit the spectral function to Euclidean solutions obtained in recent 
lattice simulations and in the gap equations formalism for gluon form factor G. The 
full dependencies of the continuous function a(u>; g(£), are indicated explicitly in 
the Ansatz H2()(l . g(£) is the renormalized QCD coupling constant at the scale £ and 
the word 'generalized' simply means that we do not assume (and do not obtain) 
the spectral function <t(uj) positive for all values of to. The spectral decomposition 
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Fig. 4. Gluon propagator for timelike momenta. 



has been fitted to some recent lattice and SDEs results (both defined in 
spacelike region) and then predicts the gluon propagator in the timelike region. For 
a details of the numerical procedure see the original paper |13| . The solutions are 
plotted in Fig.0| The spectral function has a smooth peak around p 2 = (0.7GeV) 2 
with the width « lGeV. It becomes negative for asymptotically large p 2 , as ex- 
pected already from the PT. The gluon propagator should describe the confined 
particle, so the "unusual" shape (violation of the spectral function positivity) is in 
accord with our physical expectations. 



5 Wess-Zumino model 



At the end of this paper we investigate the phenomenon of dynamical mass gen- 
eration in the supersymmetrized version of Yukawa interaction. This well-known 
Wess-Zumino model (WZM) [SUj has been already studied in the formalism of 
gap equations in couple of papers |121 140| . In both cases the authors preferred to 
work with the right number of degrees of freedom in the action and we follow this 
approach for a convenience. The auxiliary fields F, G of the Wess-Zumino super- 
multiplct (f> = (A, B, F, G, ip) were integrated out and thus we are left with the 
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following generating functional of WZM: 

Z W z[J] = J^J VADRDtj) exp / ^ ^ d 4 x [C wz + J A .A + J B .B - #] J 
iffz = Lf ree + L int (21) 
L /ree = - [d^Ad^A + d^Bd^B + ivjj flip - m 2 A 2 - m 2 B 2 - m^V] 

L mt = - mg A{A 2 + B 2 )-^(A 2 +B 2 ) 2 - gWA-igjrfripB 

where majorana fermion tp interacts through the usual Yukawa interaction with 
scalar A and pseudoscalar B. Since the whole interaction is dictated by the exact 
global supersymmctry there should not be quadratic divergences. It allows us to 
consider WZM as a natural (stable against higher corrections) low energy effective 
theory. 

Here we did not attempt to solve gap equation directly in Minkowski space, but 
apply the method first mentioned in the previous section. We will solve the gap 
equation in Euclidean space and then, by the means of spectral representation, we 
make an analytical continuation to the Minkowski space. Like in the previous study 
of gluon propagator it is based on a numerical fit of assumed spectral representation 
at Euclidean (spacelike) domain of momenta. 

In the approximation used in the paper |40j the gap equations in Euclidean 
space read: 



/•OO />7T 

n A(y) = ^2j dxx J d6sm 2 6J A 

Il B(y) = TT^ dxx d8sm 2 9J B 
27r Jo Jo 



JA = ^ ( ' GM ' n 1 2GB{X) ~ 8 x + M 2 (x) +4 (x + M 2 (x))(z+M 2 (z)) 

T(x)T{z) T(x)T(z)[y + (M(x) - M(z)) 2 } 
x + M 2 {x) (x + M 2 (x)){z + M 2 (z)) 



Jb = {SGb{x) + 2Ga{x) 



m I d9xSm9 { x + M 2 (x) IGa(z) + Gb(z)]^ 

M(y) _ 2a f°° J (« M . 2 J T(x)M{x) 

Hv) 



^j^xf^ex^e^^^-G^} (22 ) 



where x, y stands for the square of Euclidean four- momenta, z = y — 2^/xy cos 8 + x, 
and the scalar propagators are defined like Ga,b (x) = [x + Ha^x)] -1 and being 
interested in the effect of dynamical mass generation we put m = 0. Note that all 
scalar fields cubic interactions vanish in this case. For interested reader the massive 
case is reviewed and solved in |4()j . 
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Here, we instead of the use of hard momentum or Pauli-Villars cutoff we put 
higher derivative cutoff function into work. For simplicity the following form of 
cutoff function: 

dx — ► dx-7- TTTT, 7TT (23) 

(l+a:/A)(l + z/A) v ' 

has been introduced into the equations written above. The simple cutoff func- 
tion H23|) has been chosen for the purpose of convenience and we assume that any 
other reasonable choice (including Pauli-Villars regulator) gives the same results 
for p 2 « A 2 . 

There is only one scalar function characterizing scalar propagator. Thanks to 
(observed) positivity of the selfenergy II^b, the Euclidean mass function can be 
defined like Ma,b{p) = \/^a,b- The equations (|2*2*|l have been solved numerically 
by iterations. Without any doubts it provides the dynamical generation at the scalar 
sector. The dynamical masses of scalar and pseudoscalar bosons are degenerated 
which is observed even for small coupling. On the other hand, but in fact due to 
the mentioned A, B degeneracy, the majorana fermion remains exactly massless 

M(p) = 0, (24) 

which seems to be independent on the value of the coupling strength a — g 2 /4ir. 
The equation Ij24(l clearly follows from the structure of the gap equation ilL'L'l) which 
is 

M(p) ~ J dpM(q)[G A (p) - G B (p)}. (25) 

The square root of selfenergy function Ha.b is displayed in the Fig. 3a. Re- 
call, the results take sense only at the scales where p\, < A because of the cutoff 
A, further stress that the numerical results are independent on the truncation of 
momentum integrals in Il_?2[l when the appropriate hard cutoff Ah >> A. 

The observed results obviously violates the following SUSY Ward Identities : 

U A (p)=U B (p)=Tr^- (26) 

since the second equality in (|26|) is not satisfied. One can suppose that the massless 
fermion should be identified with the goldstino of broken supersymmetry. On the 
other hand we observed the similar violation also in the explicit massive case. 
To judge between the true SUSY breaking and a simple effect of gap equations 
truncation is difficult task. We could learn more from the lattice simulation, which 
hopefully will be performed in not so distinct future [4*T] . 

At this place we should mention the paper J2| where it was suggested that there 
is no DMG in the WZM. As the authors of J2| followed the same approximations 
of gap equations as here, there must some weak point in their considerations since 
there is a clear numerical evidence that II^ = Tig 7^ and here we identify this 
fault origin: The authors split the selfenergy into two pieces T1a.b(p 2 ) — P 2 [l — 
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Fab(p 2 )] + M A B (p 2 ) and they speculate that the solution M\ B (p 2 ) = should 
follows from gap equations due to the supersymmetry. Furthermore, they speculate 
that the solutions of the resulting equations for two functions Fa,b (i-e. eqs. (29) 
and (30) of mentioned paper) do not change the zero position of the propagator 
pole. However, this argument is simply wrong. Inspection by the inverse substitution 
F = ^rjjj we re P ro duce the original gap equations for II (no matter what M is, 
clearly the 'splitting' of II to M, F is meaningless here. ) 

Furthermore, one should note that scalar propagator does not posses a real pole 
but it describes resonance. Due to the behaviour of G for momenta around the 
resonant mass it is much convenient to make a fit for selfenergy function rather 
them for G directly. The most general dispersion relation reads: 

0(p) = o(A) + b(A)p 2 + [°° dx-^^. (27) 
J Q x + p 

where the function p should be non-zero from the zero threshold. Clearly fitting 
a, b, p we obtain the information about mass function M(p) in the full regime of 
momenta p. The procedure is limited by the efficiency of numerical method which 
was minimization procedure in our case. Nevertheless it offers rough numerical 
estimates when used in practice. 

Having fitted the quantities in l|27|l we have calculated the scalar propagators 
at time-like regime. In principle the matrix element of tpip — > ipip scattering can be 
read off. In the figure 3b the typical deformed Breight-Wigner shape of the square 
of propagator GG + is shown. In this case the coupling strength is a — 0.01 and the 
obtained resonant mass is m r ~ 2.3 with the width T ~ m r /4. These dimensionfull 
variables are scaled in units of 10~ 5 A 2 . 

6 Conclusion and conjectures 

We have reviewed examples of Euclidean quantum field theory models which 
provides the effect and examine the validity of standard analytical assumptions. 
Especially in the case of Yukawa theory and (strong, ladder) QED we have checked 
whether the known Euclidean solutions of Schwinger-Dyson equations are consistent 
with the assumed singularity structure of the particle propagators. Contrary to the 
soft coupling case with the explicit chiral symmetry breaking we have found that the 
obtained solutions does not posses spectral (Lehmann) representation. Although we 
do not recover the propagator singularity structure at all, the observed absence of 
propagator pole argues that such theories can not describe free particles. 

On the other hand, it is suggested that the strong coupling theory may not be 
in contradiction with standard (rather say then correct) analytical assumptions. It 
was shown that the recent lattice data on numerical simulation for gluon propagator 
can be fitted by dispersion relation with continuous non-positive spectral function. 
Here, the proposed absence of particle pole at real positive semi-axis of momenta 
is welcome and is consistent with the confinement and unobservability of free color 
objects- the gluons. 
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Fig. 5. a) Boson selfenergy b) Scalar propagators product. 



The mass generation has been studied in the WZM for the first time. Although 
the majorana fermion remains massless the nontrivial solution for the selfenergy 
functions of bosons was observed. We have also shown the evidence for resonant 
behavior of processes mediated by the scalar A and/or pseudoscalar B exchanges. 
From our numerical investigation it seems very likely that the boson propagators 
satisfy the spectral representation. 

The main results of this paper is proposal of absence of a real particle pole in 
the quantum field models with strong coupling, based on the numerical analyzes, 
the study presented here suggests that when fermion DMG effect appears then the 
'would be' particle modes condense or/and they are confined. In the certain analogy 
with QCD we can assume that only some composite objects -'bound states'- can 
freely move in a space. 

Of course, it is possible that some more realistic models of dynamical elec- 
troweak breaking (and hence mass generation of Standard model particles ) suffer 
from our simple analyzes. We have ignored particle mixing. Also the existence of 
soft and strong interacting sectors has not been properly modeled. The later prop- 
erty is a typical requirement of the Technicolors. Then, if the underlying dynamics 
is so mysterious and allows free particles to receive their masses dynamically then 
there is a white space in the road map of gap equations. Otherwise we did not 
obtain the picture which is self-consistent with the recent sophisticated nonpertur- 
bativc methods. One hopeful method of such inspection -the spectral technique- is 
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reviewed in this paper. 

References 

[1] H. Pagels, Phys.Rev.D21, 2336 (1980). 

[2] Itzykson, Zuber, Introduction to QFT; Alternatively gap equations can be read from 
the canonical formalism developed by Symanczik and reviewed in the standard text- 
book 'R.J. Rivers Path Integral Methods In Quantum Field Theory 1 . 

[3] S. F. King, Rept.Prog.Phys. 58, 263 (1995). 

[4] T. Appelquist, Lectures given at workshop 'Neutrino-Prague 2004' . 

[5] T. Appelquist, M. Piai, R. Shrock, Phys.Rev. D69, 015002 (2004). 

[6] N. Haba, N. Kitazawa, N. Okada, Invisible Technicolor, hep-ph/0504279 

[7] T. Brauner and J. Hosek, hep-ph/0505231 

[8] A. Zee, Phys. Lett. 93B (1980). 

[9] Darwin Chang and A. Zee, Phys.Rev. D61, 071303 (2000). 

[10] Kingman Cheung and Otto C.W. Kong, Phys.Rev. D61 (2000) 113012 

[11] P. Bicudo, Phys.Rev. D69, 074003 (2004). 

[12] A. Bashir, J. Lorenzo Diaz-Cruz, J.Phys. G25, 1797 (1999). 

[13] V. Sauli, submitted to Few Body Systems, hep-ph/0412188 

[14] R. Fukuda and T. Kugo , Prog.Theor.Phys. 60,565 (1978). 

[15] V. Sauli, JHEP02, 001 (2003). 

[16] V. Sauli, PhD-Thesis, available at the auhor's web pages: 
http: //gemma, ujf. cas. cz/ ~ sauli . 

[17] V. A. Miransky, Nuovo Cimento 90A, 49 (1985). 

[18] J.C.R. Bloch, PhD-Thesis, hep-ph/0208074 

[19] P.E.L. Rakow,Wuc. Phys. B 356, 27 (1990). 

[20] D.J. Gross and F. Wilczek, Phys. Pev. Lett 30, 1343 (1973). 

[21] H.D. Politzer, Phys. Pev. Lett 30, 1346 (1973). 

[22] CD. Roberts and A.G. Williams, Prog. Part. Nuc. Phys. 33 (Pergamon Press, Ox- 
ford), 1994. 

[23] R. Alkofer, L. von Smekal, Phys. Rep. 353, 281 (2001). 

[24] P. Maris, CD. Roberts, Int. J. Mod. Phys. E12, 297 (2003). 

[25] CD. Roberts, S.M. Schmidt, Prog. Part. Nucl. Phys. 45, 1 (2000). 

[26] C.S. Fischer, R. Alkofer, Phys. Rev. D67, 094020 (2003). 

[27] A. Bender, W. Detmold, CD. Roberts, A.W. Thomas, Phys. Rev. C65, 065203 
(2002). 

[28] F.J. Llanes-Estrada, P. Bicudo, Phys. Rev. D68, 094014 (2003). 

[29] J.M. Cornwall, Phys. Rev. D26, 1453 (1982). 

[30] A.C. Aguilar, A. Mihara, A. A. Natale, hep-ph/0208095 



Czech. J. Phys. -1949 (1) 



X 



17 



a to je konec tadydadydaa 



[31] A.C. Aguilar, A. Mihara, A. A. Natale, Phys.Rev. D65,054011 (2002). 

[32] A. Mihara and A. A. Natale, Phys.Lett. B482378 (2000). 

[33] J.C.R. Bloch, Few Body Syst. 33, 111 (2003). 

[34] D. Atkinson and J.C.R. Bloch, Mod. Phys. Lett. A13, 1055 (1998). 

[35] K. Kondo, Phys. Lett. B551, 324 (2003). 

[36] A.C. Aguilar, A. A. Natale, JHEP 0408 (2004) 057. 

[37] L.D. Fadeev and V.N. Popov, Phys. Lett 25B, 29 (1967). 

[38] F.D.R. Bonnet, P.O. Bowman, D.B. Leinweber, A.C. Williams, J.M. Zanotti, Phys. 
Rev. D64, 034501 (2001). 

[39] J. Wess and B. Zumino, Phys. Lett. 49B 52,(1974). 

[40] V. Sauli, submitted to Nuc. Phys. B. 

[41] M. Bonini, A. Feo, JHEP 0409 (2004) 011. 



18 



X 



Czech. J. Phys. -1949 (1) 



